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We calculate shadowing in the process of deep inelastic interactions of leptons with nuclei in the 
perturbative regime of QCD. We hnd appreciable shadowing for heavy nuclei (e.g. Pb) in the region 
of small Bjorken scaling variable 10~® < a; < 10“®. This shadowing depends weakly on Q^, but it 
may be strongly influenced, especially at a; > 10“®, by the existence of real parts of the forward 
scattering amplitudes. 
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I. INTRODUCTION 

There is considerable progress in calculations of various cross sections for deep inelastic lepton scatterii^ in the 
Lipatov limit —>■ oo and fixed at a large value Q). In particular a picture has been proposed 1 ^,^ which 
reproduces the basic features of BFKL pomeron, e.g. its Regge behavior which was earlier obtained from a direct 
resummation of the Feynman diagrams of Quantum Chromodynamics. Subsequently, the model of was applied 
to calculate the total cross section, the diffractive dissociation cross section, and the corresponding structure functions 
for the proton target [Q-H . 

In view of the successful experimental program for the so-called “small-x physics” — which is being pursued at 
HERA for proton targets and will be, possibly, extended to nuclear targets — it makes sense to adapt the results 
of [Q-|| also to nuclear targets. So, in this paper, we apply the results of [^-|| to calculate the shadowing in deep 
inelastic lepton-nucleus interactions. 

The problem of nuclear shadowing in deep inelastic lepton-nucleus interactions has been recently discussed in many 
papers. Refs. mi give an up-to-date picture of this held. 

For a description and discussion of the basic physical ideas of the model used in the present calculations we refer 
the reader to the article |^. Here a brief summary will suffice. The model was initially designed for “onium” - 
“onium” scattering, i.e. for scattering of two bound qq states with a small transverse separation. At high energy — 
in perturbative QCD — such an “onium” state evolves (cascades) into a bunch of colorless dipoles whose number 
increases as a power of the incident energy. The interaction occurs via two gluon exchange between the pairs of 
dipoles (one from each “onium”). Navelet, Peschanski, Royon and Wallon applied this model to deep inelastic 
lepton-proton scattering. They assumed — naturally — that the virtual photon at high can be described by 
an “onium”. For the target proton — which certainly does not resemble an “onium” — they made a rather bold 
assumption that it can nevertheless be approximated by a collection of “onia” with an average “onium” radius to 
be determined from the data. It was thus rather surprising that such a — seemingly crude — model described very 
reasonably the existing data on F2 and on gluon structure function in a fairly large range of and of Bjorken x. To 
obtain good fits to the experimental data the authors of Refs. |Q,|j had only to normalize F2. We show below that 
this normalization corresponds to rieff « 6.7 “onia” in a nucleon. 

In the present paper we apply these ideas to the calculations of shadowing in deep inelastic scattering from nuclear 
targets. We hope that this can provide an independent check of the hypothesis of Ref. [Q, and give more information 
about the nucleon structure at very small Bjorken x. 

The total virtual transverse photon - “onium” cross section is one of the basic inputs in computations of nuclear 
shadowing and we will use the following integral representation of = (E„P„tT„)^ (n is the index of the Good and 
Walker eigenmodes of absorption, see below) 




( 1 ) 


1 


C —200 




This integral representation is our version of discussed in |Q . We have checked through numerical evaluation of 
(H) that it agrees well with the results of Refs. y,|| provided Ueff « 6.7. In (|^) W is the number of colors involved 
in the process, and Ueff is an effective number of “onia”in the nucleon, aem = Cj is the sum of the squares of 
quark charges, and rp « .8 fm is approximately the “onium” diameter. 
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and 7 ( 7 ) is the eigenvalue of the BFKL kernel defined in terms of = ^[lnr( 7 )] as 

xh ) = 2 V'(i) - - hi ) - ■ 

x(l) determines the slope of the BFKL pomeron 
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where a is a strong coupling constant. 

Since we are interested in the structure function F 2 (crtransv + criongitud) one should add to (|^), which gives the 
transverse virtual photon - “onium” total cross section, the contribution of the longitudinal virtual photon in order 
to have a complete virtual photon - “onium” total cross section. As it was shown in Ref. Q one gets the longitudinal 
photon contribution by replacing the function h"’" {"/) of d) by h^i^y. 
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Since the integrand in ( 0 ) is strongly peaked at 7 = 1 we can approximate the above replacement by « |/i^( 7 ). 

Hence the total cross section is 
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The success of ( 0 ) (compare encourages us to apply the results of to evaluation of nuclear shadowing. 

In the next Section we present the general formula for the shadow in the deep inelastic virtual photon-nucleus 


we do the same for the “quasielastic” contributions |^]. Section |V] presents results and conclusions. Two Appendices 
give details of derivations of some formulae. 


interactions. In Section III we apply it to calculate the “triple pomeron” M contribution to shadowing. In Section IV 


II. BASIC FORMULAE FOR VIRTUAL PHOTON-NUCLEUS INTERACTIONS 


Many theoretical papers on this and related subjects have been published in the past and our task is merely to 
adapt the existing material for our purposes. Our main formula, which gives the deviation from unity of the ratio, Ra, 
of the total virtual photon-nucleus cross section to the incoherent sum of the cross sections on individual nucleons, is a 
slight generalization of the formula proposed more than 25 years ago by Gottfried and Yennie |^| and then developed 
by Yennie and Bauer et al M. It goes as follows: 
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Negative A^^ means there is shadowing, positive A^ - antishadowing. Here x = is the Bjorken variable with 
= —q^ (the negative of the four momentum squared of the virtual photon), m - the nucleon mass, v - the energy 
loss of the incident lepton (in the rest frame of the target nucleus). 
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where M is the mass of the virtual photon excitation. In this formula we neglected the correlations between the 
nucleons inside the target. The index n stands for the Good and Walker eigenmodes of absorption which 
propagate through nuclear matter without diffractive excitations, b is the impact parameter of the virtual photon, 
and Pn is the probability of realization of the n-th eigenmode. TinPn<^n = o’tot is the total transverse virtual photon- 
nucleon cross section. /„ is the forward scattering amplitude of the n-th mode. 

In |^-||] only second order interactions in dipole-dipole cross section of the diffractively excited systems were con¬ 
sidered. Then, only two classes of processes seem to be dominant Q: the “triple pomeron” contribution where the 
“onium” scatters inelastically, and the “quasielastic” contribution where the “onium” scatters quasielastically from 
the target proton. 

Note that the results of Refs. ||-||] we are going to use to calculate do not give us and /^(b) but only 
SnPr!,/rt(b)/„(b). However, S„P„/„(b)/„(b)/„(b) and higher order expressions have, so far, not been calculated 
(e.g. in the case of an exchange of three pomerons one would have to calculate a three dipoles distribution function), 
hence the last factor in (^, can, at present, be only approximated (see below). 

Note also that when the transverse spatial extension of /„ is much smaller than the nuclear radius (hence the 
transverse extension of p(b, z)) we can approximate 




(9) 


This approximation was employed in (^. Physically, this is the question of the relative transverse sizes of the cascade 
of dipoles 1^-^ and of the nucleus. We have checked accuracy of the assumption (^) for p(b, z) suitable for large 
nuclei (see Section B and found it acceptable; p(b, z)’s are taken in form of Saxon-Woods distributions. 


Now we will discuss the other input, SnPn/n/n > from 
numerous references to many formulae of the papers of Refs. 



and to make it concise we shall proceed through 


III. “TRIPLE POMERON” CONTRIBUTION TO THE SHADOWING 


We assume A ^ 1, hence a large nucleus where (||) is applicable. We first work out only the double scattering 
factor in We modify the formula (21) of ^ introducing two impact parameters b,b' (to describe collisions with 
two different nucleons) and taking the inverse Mellin transformation of both sides of (21): 
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where 


Op = [7aNcC{A) ln(l/xp)/7r] ^ = a{xp), (11) 

and denotes a partial summation over the parameters specifying diffractive eigenstates (for details see Q). The 
prime is added because the summation over the masses, M, of the diffractively excited states is left till the very end. 
W{xi 2 tXq 2 ) is a symmetric function of X 12 and a;o 2 , and e.g. for Xq 2 < X 12 


W{Xi2,Xq2) = XJ 2 
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F being the hypergeometric function, and 
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D(b, x)= f [frdz^{r, z)r In | — j exp 



where $(r, z) is the square of the proton wave function. For small Up one can evaluate ( 0 ) analytically (the relevant 
integrals are given in Appendix A) and obtain for the transverse virtual photon 
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where 


yil) = J llA-^'l',P,y‘^)dy , (15) 
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and ro = f d^rdz‘l>(r, z)r is a non-perturbative parameter determined from the data. A fit to the proton structure 
function gives rg ~ .8 fm 

At this point we introduce the nuclear algorithm with A collisions inside the nucleus. We write 
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The question now is what to take for f^ff ? As we have already remarked the results given in do not go 

beyond the double scattering. Also, do not provide us with amplitudes for the eigenstates of absorption but 

with amplitudes averaged over many eigenstates. In order to have such amplitudes one would have to do new and 
— possibly — rather complicated calculations. Therefore, in order to estimate the multiple scattering corrections we 
shall use the full amplitude for dipole-proton scattering at fixed p (this is the diameter of the “onium” representing 
the virtual photon) 
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The integration over d^b can be done (see Appendix A) and we obtain 


feff{ro,p,Xp) = 2^/2 
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A comment is here in order. To estimate the multiple scattering corrections one can use some other than feff 
effective amplitude, e.g. the one proposed in Ref. Eq. (5). In our notation it amounts to taking the ratio 
{fn)lfeff instead of Jeff- As expected, (/^)//e// > feff but both are of the same order of magnitude. Since the 
complete expressions for multiple scattering corrections need direct calculations of the contributions of interactions 
with 3 and more nucleons which are not done, we choose feff as a somewhat simpler expression. 

The last point is to average Ta over the probability, $, to find a dipole in the virtual transverse or longitudinal 
photon Q: 
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Here 


Q = ^/z(l - z)Q 


( 22 ) 


Thus the “triple pomeron” (TP) transverse (T) and longitudinal (L) contributions to the shadow, within the ap¬ 
proximation (^, is as follows 
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For comments on the limits of the integration over (3 and p see Section ^ 

It is convenient to extract from ( p^ a correction factor, Cmsix,Q'^), for multiple interactions higher than double. 
This is the factor which multiplies the double scattering contribution to give A^'^(a;, Q^). Thus we have 
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Crasix, Q^) will be used in the next section to compute the shadow of the whole multiple scattering process from the 
shadow of the double scattering contribution. We believe that the level of accuracy of our calculations and the size 
of Cmsix,Q^) justify this approximation. 


IV. THE QUASIELASTIC CONTRIBUTIONS TO THE SHADOWING 

This time the Ref. 0 provides us with the basic formulae and we will assume A 1 hence, again, we accept the 
approximation i)- The method used in ||^ to obtain the double scattering contribution is different than in the case 
of “triple pomeron” and we now have at our disposal the probability amplitude for the intermediate state which, 
in the present case, is the one dipole state. Therefore, the procedure of integration over b is done at the level of the 
amplitude which, subsequently, is contracted with the probability amplitude for finding a dipole in a photon of a given 
polarization, and the result squared. 

So, starting from the formulae (19) and (20) of Ref. @], we hrst replace '^^^2 by ^ integration 

of the amplitude T over b, and obtain the following double scattering terms in the approximation (|^): 
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where feff{ro, P, Xp) is given by (|^, and M = y/IiT 


z)M, with M being the mass of the virtual photon excitation. 
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In the double scattering cross sections ( pS]) Jeff appears quadratically. Therefore, since the basic ingredient of 
Cms(x, Q^) defined in Section HI is the factor [1 — /e// dz^pib, 23 )]"^“^ (defined through the same amplitude Jeff 

as in dll)), we accept that by multiplying for just the double interaction inside the target nucleus by Cms we 
obtain a reasonable approximation for with all possible (i.e. A) interactions included. The error made in this 
approximate procedure is not significant because Cms is close to unity throughout the relevant region of x and 
(see Fig. 3). 


Following steps similar to those which led us to the formula (|23D we obtain for the quasielastic contributions to the 
shadow of the transverse and longitudinal virtual photons the following expressions: 
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where feffYo.p.Xp) is given by 

Clearly, setting Cms = 1 reduces both (U) and (^) to the correct expressions for the double scattering contributions 
to the shadow. In fact — as we have already stressed — the double scattering term is more reliable than the multiple 
scattering correction because its input from j6|-||] comes exclusively from the “first principles”. That is why it is of 
interest to discuss deuterium target. This will be discussed in a forthcoming paper. 


V. DISCUSSION OF THE RESULTS AND THE CONCLUSIONS 

Before presenting any results a comment on the real parts of the amplitudes / is in order. First, the dependence 
of atot on X , hence on 2mi' « s implies that Re[/] must be different from zero. This follows from the derivative 
analyticity relations for the forward elastic hadronic amplitudes worked out in e.g. Ref. [Q. We can introduce, 
approximately, this effect into all our expressions derived so far performing the following replacements: 

(/5^(/n)(l-*a)^ 

feff /e//(l - *d), (27) 

where a is the ratio of the real to imaginary parts of the forward amplitude. We use the leading term of the expansion 
of the forward amplitudes / given in [Q : 


a = 


Re/ 

^tot 


tan (^TrAp) 


(28) 


We accept Ap .3, hence d « .5 . 


Now the results. We have calculated all contributions to the shadow A/i from the triple pomeron and the quasielastic 
diffractions for two target nuclei: Lead (Pb,A=208) and Calcium (Ca,A=40). All integrations in (^) and (p^ ) were 
done numerically with an extensive help of Matheinatica (more details are given in Appendix B). 


The single nucleon densities p(s, z) were taken in the standard form 
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with r = \/s^ + z"^, ta = 1.2 and a = .54 fm. 


V = 




(30) 


so that f d^sdzp(s, z) = 1. 

We assumed Ap = .3 which , through relation (^) and using = 3 gives the strong coupling constant a = .11. 
This is an acceptable set of parameters characterizing the BFKL pomeron. We set Ueff = 6.7 to recover the fit to 
the experimental data of the F 2 of the nucleon [^|-||. 

In all integrations we restricted the values of Xp to 


= ^<.1 




(31a) 


which, as seen from (^, amounts to excluding distances {zi — Z 2 ) < 2fm from the double- and higher multiplicity 
interactions; physically a reasonable condition to impose. This was done because in construction of the triple pomeron 
contribution to A^i, Eq. (^), see also Appendix A, we expli citly used the approximation of small a 
We have checked that changing the upper limit for Xp 


,p, hence small Xp. 
(|31a|) by a factor 1.5 or 2 does not change significantly 


at the values of x for which we have done the calculations, namely x < .01. This gives the lower limit for integration 
over (3: P > IQx. 

The upper limit for fd on the other hand, should be taken close to unity because small masses give important 
contributions to shadowing. So, the integration range for [3 is taken as 


lQx<P<l . 


(31b) 


All this makes almost all asymptotic formulae for the “triple pomeron” given in Refs. ||^-^ inacceptable for our 
purposes because they lead to infinities for integrations over /3’s close to unity . 

We wish to make it clear: by extending our calculation to P close to 1, hence to very small masses M, we go beyond 
the region of /3’s for which the original formulation of Refs. was devised. We make this extrapolation because 
the complete integral representations of various TinPnfnfngiven also in |6|-|^ lead to well defined finite results, even 
for (3—^1 and only they are employed in our calculations. The price to pay is loss of simple analytic expressions 
and necessity of numerical integrations. Although the asymptotic form of atot = ^nPnO'n does not lead to any such 
problems, to be consistent with our treatment of , we also used its integral representation (^. 

We set the cutoff parameter, r*, of the integration over p, following the discussion of Ref. 0, at r* = tq. Unfor¬ 
tunately, in the case of nuclear targets the saturation of the values of the integrals over p does not look as good as 
in the case of the nucleon target (see Ref. 0, Fig. 2). Therefore, we should treat r* as a parameter with the lower 
bound equal tq. This implies that taking r* = vq we tend to underestimate the shadow. 

One more general comment: we accept that the contribution to A^ of the triple pomeron - , of the quasielastic 
transverse - , and of the quasielastic longitudinal processes are additive. In other words we do not introduce any 
interferences between them. 


In Fig. 1 a,b the total shadow A^ = -|- + A^'®^ -1- is presented for Pb and Ca targets. As we 

can see, the pattern of all curves for these two targets is very similar, only the size of the shadow for Pb is a factor of 
two larger than for Ca. We have also checked through explicit calculations that all components of A a of Pb are by 
the same factor (« 2) larger than of Ca. Therefore, we will, henceforth, show only the results for Pb. 


In Fig. 2 a,b we show A^^ with (a) and without (b) presence of the real part of the forward scattering amplitudes. 
Comparing this figure with Fig. 1 we see that A^^ dominates A^, and that the presence of d 0 makes shadowing 
smaller. In fact a different from zero real part makes |Aa| smaller: by 16% at very small x and by 60% (!) at a; = .01 
. (The role of the real parts of the forward amplitudes at various x has been recently discussed in 14|l^.) 


In Fig. 3 a,b we show the correction factor for multiple scattering higher than double, Cms{x,Q^), with (a) and 
without (b) presence of a. We can see that in either case C'ms is close enough to unity to justify our procedures in 
construction of A_ 4 . 


7 







1 2 3 4 5 



ln(Q^/GeV^) 


FIG. 1. The total shadow, Aa, for (a) Pb and (b) Ca nuclei. —Aa is plotted vs ln(Q^/GeV^) 
scaling variable x. The ratio of the real to imaginary parts of the forward amplitudes is set at d 


for four values of the Bjorken 
= tan(i7rAp) « .5 . 


The following conclusions result from inspection of the Figs. 1 — 3: 

(i) In the region x < 10“^ our calculation shows an appreciable shadowing, increasing with decreasing x. The x 
dependence of Aa follows approximately the power law 

Aa ~ (32) 

with A « 0.25, somewhat smaller than the assumed value of Ap = 0.3. An indication of saturation is seen at 
q 2 =4GeV^ 

(ii) A A decreases slowly with increasing Q^. The effect is somewhat stronger at small x and large Q^. 

(iii) Multiple scattering corrections account for less than 20 percent of the value of Aa- 

(iv) The existence of the real parts of the forward amplitudes makes shadowing appreciably smaller, especially at 
X close to X = 10“^. 
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FIG. 2. The “triple pomeron” contribution to A a with a 7 ^ 0 (a), and with d = 0 (b), for Pb target. 


Several comments are in order: 

(a) The rather weak dependence of Aa seems puzzling if one takes into account that in the dipole model 
the total (virtual) photon-hadron cross section decreases as Q~^ at large Q |^,^,||. One would thus naively expect 
Aa ~ (/^)/(/) to decrease also as Q~^ which is another version of the Bjorken-Gribov @0 paradox. In the dipole 
model the paradox is resolved (and approximate scaling of Aa recovered) due to the cascade nature of the intermediate 
state which interacts inside the nucleus. This cascade nature implies very strong correlations between the dipoles in 
the onium representing the virtual photon so that the double dipole density does not resemble at all the product of 
single densities and both (/^) and (/) behave as Q~^ (apart from logarithmic corrections). The approximate scaling 
of Aa follows. It is interesting to observe that this explanation is radically different from that proposed in the parton 
model (first paper of i, 

(b) It should be pointed out that, although our calculation is based on perturbative QCD (as is the whole dipole 
approach), it nevertheless requires two “non-perturbative” parameters (rp and Ueff) describing structure of the 
nucleon. These parameters were estimated from the data on the total (virtual) photon-proton cross section as they 
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cannot be calculated in the dipole model. Particularly important for the absolute value of the shadowing turns out 
to be rieff, the average number of “onia” forming the proton. This confirms the point of view that the perturbative 
QCD alone cannot describe fully deep inelastic scattering even at very small x 1^. It is rather remarkable, however, 
that the required “non-perturbative” input is so restricted. 

(c) As was already pointed out in the main text, our calculation relies heavily on the calculation of the deep inelastic 
diffraction reported in . This is a consequence of close relation between the nuclear shadowing and diffraction (see 
e.g. Ip^ ). Nevertheless, it is worth to mention that — due to the power-law tail in the impact parameter dependence 
of the quasielastic and triple pomeron amplitudes — nuclear shadowing is much more sensitive to the details of the 
effective cut-off in the photon-nucleon impact parameter plane than is the total diffractive cross section. At this point 
the two calculations are substantially different. 

(d) Our treatment of multiple scattering is only approximate. To improve on this point, it is necesarry to calculate 
directly at least the contribution of interaction with 3 nucleons which would give an estimate of the accuracy of our 
approximation. This seems a feasible, although probably a complicated calculation. 
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FIG. 3. The correction factor Cms{x,Q^) for higher than double scatterings with d 7 ^ 0 (a), and with 0 = 0 (b), for Pb 
target. 
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APPENDIX A: 


Here we give more details of some intermediate steps of calculations of Section. 3. First the integration over b and 
in It can be done using the formula 


ln(—)e = TT +/i)e = —e 


(Al) 


max(r,K) 


where h = ln(r/a:). 

The integration over cPrdz is done with the help of the (asymptotic, Op —> 0 ) formula 
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where rg = f ^(r)rd^rdz. The argument justifying (K^) goes as follows. We start with the identity 
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[ $(r)re ^^'^d'^rdz = 2.[-^Xp^x [ 

J y 2ap P J 2TTi 


d'J aNc 


Xil) ln(a:p) 




d^rdz 


Let us denote 


J ^{r)r^d'^rdz = 


(A3) 


(A4) 


(A5) 


Of course fg depends on 7 . But, in the limit Op —*■ 0, the main contribution to t he integral (A4) comes from the 
region 7 « 1 and we approximate fg by rg = ©(7 = 1). Using this we obtain from (A4) 


J y 2ap J 2Tri ^ x 


m 


We use again (^) and obtain ([A^). 

Now the formula &■ Introducing ([Al|) and (|A2| ) into (^^ we obtain 

/ T~ n~, , - ~ l6TratffaNcrg /a.\ 

dHdH\KPnfnih)fnm) = -^ (^) 

Xp 


C+200 


_ p2-^ p-aNax(l) /tt 


2TTi 


X J dxi2dxg2W{xi2,X02)e 2“p('*0+^l) 


(A7) 
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where ho = In(^), hi = lii(^) . 

Next point is the integration over dxo 2 dxi 2 . This we do as in 1^ (although the accuracy of this procedure may be 
not very good). Let us repeat the steps. Since the integrand is symmetric in a:o 2 , a;i 2 we have 


OO OO Xi2 

1= [ dxi2dxo2fL(a;i2,a;o2)e"5“p('*o+^?) = 2 [ dxi2e~^'"^’"^^xj^^ [ — ) . (A8) 

J J J Xi2 Xi2 


The second integral can be approximated (in the limit Up « 0) by e 2°'p^^V('y), where 1 ^( 7 ) is given below Eg. (p^ . 
Substituting this into (A8) we obtain 


OO 

J Xi2 


(A9) 


This is an easy integral (gaussian, after substituting u = ln(^)) and we obtain finally 


/ = 2E(7)rJ 



This completes the integrations in (A7) and we obtain Eg. 




(AlO) 


APPENDIX B: 

We comment here on our numerical calculations. The contour integrals in Egs. (P and (|lj) and the integral of the 
hypergeometric function defining V ( 7 ) were calculated with Mathematica. 

For a fixed value of Bjorken's x, the contour integral in © becomes a function of two variables: p and /3. Such 
a two-dimensional function was tabulated (for x = = 10 “"^, x = 10 “^ and x = 10 “^) and used as an input in 

the Fortran program calculating the triple pomeron contribution (|2^). Similarly, the integral over z of the probability 
distribution $(£, p, Q') was first tabulated with the help of Mathematica, and then used as an input in the Fortran 
program calculating (^). 

Mathematica was used once again to calculate and tabulate the integrals over p of the Bessel functions in Eg. (|^). 


[1] V.S. Fadin, E.A. Kuraev and L.N. Lipatov; Phys. Lett. B60, 50 (1975); L.N. Lipatov, Sov. J. Nucl. Phys. 23, 642 (1976). 
For a recent review and further references see L.N. Lipatov, Acta. Phys. Pol. B26, 1815 (1995). 

[2] N.N. Nikolaev and B.G. Zakharov, Z. Phys. C49, 607 (1991); C53, 331 (1992); C64, 651 (1994); J. Exp. Theor. Phys. 
78, 598 (1994). 

[3] A.H. Mueller, Nucl.Phys. B415, 373 (1994); A.H. Mueller and B. Patel, Nucl. Phys. B425, 471 (1994); A.H. Mueller, 
Nucl. Phys. B437, 107 (1995). 

[4] H. Navelet, R. Peschanski and Ch. Royon, Phys. Lett. B366, 329 (1995). 

[5] H. Navelet, R. Peschanski, Ch. Royon and S. Wallon, Phys. Lett. B385, 375 (1996). 

[6] A. Bialas and R. Peschanski, Phys. Lett. B378, 302 (1996). 

[7] A. Bialas and R. Peschanski, Phys. Lett. B387, 405 (1996). 

[8] A. Bialas, Acta. Phys. Pol. B27, 1263 (1996). 

[9] General reference: “Light and Heavy Nuclei in HERA”, in “Future Physics at HERA”, Proc. of the Workshop 1995/96, 
Vol. 2, ed. by G.Ingelman, A. De Roeck and R. Kalmer, DESY Sept. 96. 

[10] N.N. Nikolaev, B.G. Zakharov and V.R. Zoller, Z. Phys. A351, 435 (1995). 

[11] G. Piller, W. Ratzka and W. Weise, Z. Phys. A352, 427 (1995). 

[12] N. Armesto and M.A. Braun, “The EMC effect at low x in perturbative QCD”, preprint US-FT / 15-96, eprint hep- 


ph/9603360. 


12 









[13] 


M. Arneodo, A. Bialas, M.W. Krasny, T. Sloan and M. Strikman, “Nuclear Beams at HERA”, in “Future Physics at 
HERA ”, Proc. of the W orkshop 1995/96, Vol. 2 p. 887, ed. by G.Ingelman, A. De Roeck and R. Kalmer, DESY Sept. 96, 


[14] 

[15] 

[16] 


[17] 

[18] 

[19] 

[ 20 ] 
[ 21 ] 

[ 22 ] 

[23] 


eprint hep-pli/9610423 


M. Strikman, Acta Phys. Pol. B27, 3431 (1996). 

K. Gottfried and D.R. Yennie, Phys. Rev. 182, 1595 (1969). 

D.R. Yennie; Acta. Phys. Pol. B7, 897 (1976); T.H. Bauer et al. Rev. Mod. Phys. 50, 261 (1978); see also: G. Grammer 
Jr and J.D. Sullivan, in “Electromagnetic Interactions of Hadrons”, ed. by A. Donnachie and G. Shaw, Plenum, New York 
(1978) Vol. 2. These review articles contain references to the earlier papers, notably to the seminal papers by L. Stodolsky. 
M.L. Good and W.D. Walker, Phys. Rev. 120, 1855 (1960). 

D.P. Sidhu and U.P. Sukhatme, Phys. Rev. Dll, 1351 (1975). 

A. Bialas and W. Czyz, Phys. Lett. B328, 172 (1994). 

V.N. Gribov, Sov. Phys. JETP 30, 709 (1970). 

J.D. Bjorken, Lecture Notes in Physics 56, “Current Induced Reactions”, ed. J. Koerner, G. Kramer and D. Schildknecht, 
Springer-Verlag (Berlin, 1976), p. 93; J.D. Bjorken and J. Kogut, Phys. Rev. D8, 1341 (1973). 

L. L. Frankfurt and M. Strikman, Nucl. Phys. B316, 340 (1989). 

J.D. Bjorken, “Rapidity Gaps in Deep Inealstic Scattering”, preprint SLAC-PUB-7096 (1996), eprint hep-ph/9601363 


13 




